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Abstract: We give a lattice theory treatment of certain one and two dimensional quantum 
field theories. In one dimension we construct a combinatorial version of a non-trivial field 
theory on the circle which is of some independent interest in itself while in two dimensions 
we consider a field theory on a toroidal triangular lattice. 

We take a continuous spin Gaussian model on a toroidal triangular lattice with periods 
Lq and Li where the spins carry a representation of the fundamental group of the torus 
labeled by phases u$ and u\. We compute the exact finite size and lattice corrections, to the 
partition function Z, for arbitrary mass m and phases Ui. Summing Z -1 / 2 over a specified 
set of phases gives the corresponding result for the Ising model on a torus. An interesting 
property of the model is that the limits m — > and Ui — > do not commute. Also when 
m = the model exhibits a vortex critical phase when at least one of the Ui is non-zero. 
In the continuum or scaling limit, for arbitrary m, the finite size corrections to — InZ are 
modular invariant and for the critical phase are given by elliptic theta functions. In the 
cylinder limit L\ — > oo the "cylinder charge" c(mo,w 2 Lq) is a non-monotonic function of 
m that ranges from 2(1 + 6uq(uq — 1)) for m = to zero for m —>■ oo but from which one 
can determine the central charge c. The study of the continuum limit of these field theories 
provides a kind of quantum theoretic analog of the link between certain combinatorial and 
analytic topological quantities 
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§ 1. Introduction 

In this paper we examine some discretely formulated quantum field theories as well as 
successfully computing their continuum limits. We shall consider field theories in one di- 
mension, over the circle S , and in two dimensions, over the torus T 2 . However a crucial 
feature of the fields will be that they will not be singly-periodic and doubly-periodic func- 
tions over S 1 and T 2 respectively; instead they will be sections of appropriate bundles over 
S 1 and T 2 . These sections will then not be periodic but will possess a non-trivial holonomy 
when rotated round a non-trivial loop on their base space. A brief summary of some of 
these results has been presented in [1] . 

In our two dimensional model we shall be able to carry out an extensive calculation 
of finite size effects which are associated with the onset of continuous phase transition. 
These effects provide corrections to the bulk system behaviour and play a vital role in 
our theoretical and experimental knowledge of systems at and near criticality cf. [2]. 
Such effects have a complex interplay with the other features of these systems giving rise 
to phenomena such as crossovers from one characteristic behaviour to another [3]. Two 
dimensional physics has provided a framework where substantial progress has been made 
in understanding such problems and considerable interest has been focused on them. Some 
of these calculations use generalised Gaussian models [4,5] and our work can be considered 
in a similar spirit. 

In one dimension we use a zeta function to regularise the theory and also construct 
the theory discretely on a polygon. We find that when the Van Vleck Morette determinant 
is included in the measure of the discretely calculated partition function the massless case 
is independent of the number N of points in the polygon. In addition we find that the 
combinatorial and continuum partition functions agree precisely and are equal to a power 
of the Ray-Singer torsion for the appropriate flat bundle of sections over S . 

For our two dimensional theory we obtain the exact finite size corrections to the free 
energy of our lattice model. These corrections form a modular invariant expression and 
for the critical phase this expression is expressible, in the continuum limit, in terms of 
elliptic theta functions and the associated partition function holomorphically factorises in 
the modular parameter r. The model possesses two holonomy phases (cf. section 3) u$ 
and u\ and a mass m; there is a critical vortex phase when the mass is zero and one of the 
phases is non-zero. Modular invariance persists whether m is zero or not. 

The principal new results contained in this work are 

(i) A demonstration that the combinatorial torsion arises naturally when the Van Vleck- 
Morette determinant is incorporated in the integration measure. 

(ii) The determination of the exact partition function for a massive Gaussian model where 
the spin variables are sections over a generic triangulated torus. 

(iii) The demonstration that the finite size corrections are modular invariant irrespective 
of the value of the mass and the demonstration that in the massless case they are 
related to a topological quantity known as the <9-torsion. 

(iv) The determination of equivalent data for the two dimensional Ising model on a trinan- 
gular lattice with three independent nearest neighbour couplings. 

The organisation of the material in the rest of this paper is as follows. In section two 
we discuss our one dimensional theory in its continuum and in its discrete formulation. 
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Section three moves on to two dimensions where we formulate a gauge theory based on flat 
bundles over the torus T 2 and then construct the associated lattice model. This section 
contains some of our basic results on modular invariance both at and near the critical 
phase. In section four we deal with bulk properties the lattice model including those due 
to finite size and finite lattice spacing. We also examine the conformal properties of the 
model. Finally section five contains our concluding remarks. 

§ 2. A twisted lattice theory in one dimension 

We start in one dimension with a gauge theory on a circle S 1 of circumference L. Any 
connection A of such a theory necessarily has zero curvature; nevertheless A can still be 
non-trivial since it can have non-trivial holonomy. The action S of our theory is given by 

S=\f (fd* E d E f + m 2 ff) (2.1) 

where (as mentioned in the introduction) / is not a function on 5" 1 but a section of a flat 
bundle E over S 1 , finally oIe is the covariant derivative which acts on such sections. 

A few details about the flat bundle E: The sections of E must have non-trivial holon- 
omy and, for simplicity in the subsequent discretisation of the theory, we take the holonomy 
group to be U(l). Now a general, rank n, vector bundle E over S 1 is a sum of appropriate 
powers of line bundles Ci i.e. 

E = Cf 1 ®C* 2 ®...©£°n j witha.GZ (2.2) 

where Ci is a single line bundle over S 1 . 

Because a completely general E is expressible in terms of line bundles Ci we shall just 
treat the special case where E is a single line bundle and take E = C a section / of which 
has the property that 

/(x + nL)=e 2 ™7(x). (2.3) 

The computation of the partition function for a general E (of rank n) is an easy extension 
of the line bundle case. 

Now the partition function Z(S 1 , £, m) for the action S above is given by the functional 
integral 



Z(S 1 ,C,m) = JvfVfexp -± J f(d* c d c +m 2 )f 



(2.4) 



This is immediately calculable (cf. [6] and [7] and references therein) and is given by 

Z(^,Am) = |det[ ^+ m2 l| 1 (2.5) 

where \i is a mass scale introduced to render the partition function dimensionless. In 
the above equation &et[ d * cdc J[ rn ] is given a meaning by defining it through the (- function 
regularization procedure. To implement this procedure we need to find the eigenvalues of 
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the operator d* c dc + m 2 . Eigensections replace eigenf unctions and a basis of eigensections 
is given by 

e n (z) = -^exp 2 ^^ (2.6) 
Now if we define the inner product between sections by 

<f,g>= f L dxf(x)g(x) (2.7) 
Jo 

The relation < dcf,g >=< /, d* c g > defines the action of d* c and we find the eigenvalues 
are 

An = ( Mn + «) )2 + m2 (2g) 
The ^-function of interest to us is then 

Cd*d c +mi{s) = (^— ) 77 — ; — 777~7777I7777"- ( 2 - 9 ) 

The easiest way to evaluate the sum is to use the Plana sum formula [8] which can be cast 
in the form 

v i = -* 1/2 ns-\) _ s+ i 

^ ((n + u) 2 + x) s r(s) X 

sin(Trs) f°° J ___ 2a d 



H — - / dqq Zs — In 

7T J dg 



■y g — 1~K\J q 2 +x+2iriu 



(2.10) 

From which we deduce that 

Cd*d £ w(0) = (2.11) 

and 

C^ d£+ma (°) = - mL - ln I 1 - e" mL+27riu | 2 . (2.12) 

Notice that the arbitrary undetermined scale \i has completely disappeared from the answer 
and we are led to the surprisingly elegant result that the 

and hence T = — \n.Z{S 1 , £, m) is given by 

V = mL + In \l-e- mL+27Ttu \ 2 (2.14) 

In the massless limit, m = 0, we find that Z(S 1 , £, 0) is independent of the circumference L 
of S 1 , a reflection of metric independence. In fact the resulting expression is a topological 
invariant of the bundle 

..... 1 . 
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where T(S 1 , L) denotes the Ray-Singer analytic torsion^ of the line bundle C and is given 
by 

T(S\C) = {2sin(7ra)} 2 . 

Though the resulting expressions above are elegant many aspects of the procedure that 
leads to them are mysterious. For example, from general principles one would expect that 
the partition function would have a term proportional to the number of degrees of freedom, 
which would give rise to a divergent contribution in the continuum limit. It is therefore 
illuminating to examine the above theory from a lattice point of view. We therefore turn 
to its discretisation. 

To this end we replace the circle 5" 1 with its natural combinatorial counterpart namely 
a polygon with N vertices. The discretised basis of eigensections is then take to be 

{— L e * x " fc : k = 1, . . . , N} where x n = ^ + ^ , n = 0, . . . , N - 1. (2.16) 
V Na N 

where the lattice spacing a = ^. Since the connection is locally trivial the action of the 
discretised derivative dc on sections can be taken to be 

dcm = MJS (2 . 17) 

ct 

The discretised inner product becomes 

(JV-l) 

<f,g>=a f(k)g(k) (2.18) 

k=0 

and hence 

d*/( fc ) = (2.19) 
a 

Applying these rules to the basis e n (k) implies that 



dce n (k) = -^e 1 ? sin e n {k) 
and d* c e n (k) = ——e~ 1 ^ sin e n (k) 



(2.20) 



and so we find that the finite dimensional lattice version of det[d^dc + m 2 ], which we shall 
denote by det jy[d* c dc + m 2 ], is given by 

detN[d*£dc + m 2 ] = Y\ A n , where A n = — sin(x n ) 2 + m 2 (2.21) 



a' 

n 



(°) xhe equality of the combinatorial and the analytic torsion for C provides one with a topological 
derivation of Lerch's classical formula: Q' U {Q) = —2 In 4 — 4 In sin(7ru), where Cu(s) = y^^__^2/ (n + u) 2s . 
In fact this is pointed out in Cheeger's paper [9] though the details are different: Cheeger considers, not a 
(complex) line bundle C, but a (real) SO(2) bundle E. 
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The corresponding lattice partition function we denote by Zn(S 1 , £, m), and it is, of course, 
automatically finite. We have 

Z N {S 1 ,C) = I Jjd0i(fc)d0 2 (fc) expf-S] where S = ^ < d c p,d c p > . (2.22) 

A; 

We have taken = 0i + z</>2 and </? = \/a0, then </? corresponds to the usual continuum 
field and is a dimensionless lattice field. This identification renders the partition function 
dimensionless as we require it to be and in fact we have 

Z N (S\C, m) = {^ ( ^+ m ' )aa }"' (2.23) 

Defining W = — In Z N we find that W is given by the sum 

A„.a 2 



W = ln 



2tt 

n=0 



N-l 9 9 



(2.24) 



- AHnTT + ^ ln[l H cos(x n )] 



n=0 

The sun can be performed using the basic identity 

JV-l 



ln[z - cos(a; n )] = iV ln[^±] + ln |1 - z V™T (2.25) 

n=0 

where z± = z ± V^ 2 — 1- Hence we find that W splits in the form 

W = NW B + W F 

where the extensive term is known as the bulk contribution. The quantity Wb is the 
contribution to W per lattice site in the thermodynamic limit i.e. 

W 

W B = hm — 

and is given by 



m 2 a 2 m 2 a 2 + rrficfi 



W B = ln[l + — ^— + y - ] - In 2tt. (2.26) 



The remaining term Wf captures the effects of a finite lattice and is given by 



W F = ln |1 - (1 + — - Jm?a? + ^)^V™f (2.27) 

Z V 4 
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The continuum or scaling limit is a constrained thermodynamic limit and corresponds 
to taking N — > oo while keeping Na = L fixed. In this limit we have 

NW B = -Nln2n + F B where T B =mL (2.28) 

and 

r F = ln|l-e- mL+2 ^ m | 2 (2.29) 

We see therefore that the quantity W + N In 2n tends to a finite continuum limit which 
in fact coincides with the result provided by the (^-function prescription for the functional 
determinant. 

It is clear from the above that we could alter the integration measure to absorb the term 
— iVln27r. This is most naturally done by including the Van Vleck-Morette determinant J 
in the integration measure so that 



Z c N {S\C,m) = / n#i(/c)# 2 (/c)Jexp[-S] 



where J = < det 



d 2 S 



2n dfaifydfaik - 1) 
The resulting partition function Z^ f (S 1 , £, to) is given by 



1/2 



(2.30) 



Z c N (S\C,m) = z%\l-z^e 27ri 

, m 2 a 2 I : : m% 4 ( 2 - 31 ^ 

where z± = 1 H ± \ m z a z H — 



With L = Na fixed we have 



lim Z C N (S\C, to) = Z(S\ C, to) (2.32) 
But notice that if to = then 

Z »( Sl > C >°)=T(&-Zj (2 - 33) 

In other words the discrete partition function with the Van Vleck-Morette determinant 
included in the continuum limit yields exactly the expression obtained by the ^-function 
method; in addition in the massless limit this lattice partition function Zj^(S 1 , £, 0) coin- 
cide exactly with Z(S 1 , £, 0) and is therefore inversely proportional to the torsion without 
the need for a limit. The partition function Z^{S 2 , £, 0) is therefore closely related to the 
combinatorial torsion and can be viewed as providing a realization of this object. 

§ 3. A two dimensional theory at genus one 

Now we go to two dimensions where we want to study a field theory on a Riemann surface 
of genus one i.e. a torus which we denote by T 2 . We shall then go on to discretise with a 
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specific triangulation to be described below. We begin by we providing a brief summary 
of the continuum setting. 

First of all let us realise the continuum torus T 2 as a complex manifold using the 
standard quotient 

T 2 = c/r (3.1) 

with T = Z © Z ~ tti(T 2 ) (3.2) 

Next, because we are going to have flat connections, we need a flat bundle E over T 2 , this 
is obtained by choosing an element x °f 

Hom(7r 1 (T 2 ),G) (3.3) 

where G is the bundle structure group. We take E to be a line bundle C so that G = U(l). 
Then X e Horn (tti(T 2 ), G) is of the form 

x:ni{T ^^S 1 (3.4) 

and x is then actually a character of iii(T 2 ) which is the reason for adopting the present 
notation. 

The action we choose is given by 

S=\f Hd* c d c +m 2 )is (3.5) 

where m is a mass. Now the torus is a Kahler manifold and so has the property that 
d* c dc = 2d* c dc, where d = d/dz. This means that the partition function Z(T 2 ,£,m) for 
this action gives the <9-torsion if we set m = 0, i.e. 



Z(T 2 ,£,0) = Jvi/Di/exp ~J vdcdcy 



d* c d c ,\ (3.6) 



det(^) 
= {T S (T 2 ,£)}- 1 

where Tg(T 2 ,C) denotes the <9-torsion( b ) and det(d^dc) is calculated (as is the analytic 
9-torsion) using a zeta function and it is useful to note in the present context that 

d* £ d c = 2d* c B c 
=► Ckdd ) = Cka £ (°) " ln(2)C^a £( o)(0) 

(3-7) 

= C^*a £ (0), since calculation will show that C^a^o)^) = 
det(d* c d c ) = det(d* c d c ) 



^ The 9-torsion is the complex analogue of the ordinary Ray-Singer analytic torsion T(M, E). For even 
dimensional M the ordinary torsion is always unity; however if the manifold M is complex then there exists 
the analytic 9-torsion which is non-trivial cf. [10]. 
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Now the d-torsion T a (T 2 , C) is a rather complicated object and was shown in [10] to 
be expressible in terms of a ratio containing a Jacobi theta function and the Dedekind eta 
function. Before giving Tg(T 2 , C) we need to give a parametrisation of the homomorphism 
X and to do this we use a notation similar, but not identical, to that of [10] : Let the torus 
T 2 be denoted in standard fashion by the usual parameter r G C + (C + denotes the upper 
complex half-plane Imz > 0), the generators of wi(T 2 ) are denoted by the pair {1, r} and 
an element of 7Ti(T 2 ) by the pair of integers (m, n). Then if u , u\ G [0, 1] one has 



(m, n) h 

While for the T^(T 2 , C) we have 



^S 1 

exp[27ri(muo + nu\)} 



(3.8) 



. 2 ,0i(iti - r« ,r) 
exp|7rm 7~J — 



rj(r) 



(3.9) 



We finish this continuum summary with a description of the continuum spectral data 
for d^dc'- For simplicity we shall work with the Laplacian d*^dc and denote it by A^; 
its eigensections by £ mn (x,y) and eigenvalues by A mn can be computed fairly straightfor- 
wardly. The details are as follows: let z = x + iy denote the local coordinate in C then 
take the torus to have periods L and L\ but at an angle 6 to one another. Then a point 
z and z + (mr + n)Lo represents points that are identified under the quotient action (3.1) 
where 

with Tn = ^-cos6> and Ti = ^sin#. (3.10) 



r = r + ITi 



It is sensible to change (x,y) to new (real) variables (xq,xi) defined by writing 

z = xq + x\ cos 6 + %x\ sin 9. 

In these new co-ordinates the metric ds 2 = dx 2 + dy 2 becomes ds 2 = dx$ + 2 cos Odxgdxi + 
dx\ and the Laplacian Ajr becomes 



sin 2 6 



d 2 



Oxq 



— cos 9 



d 2 



+ 



d 2 



dxodxi dx\ 



It is then easy to find that the eigensections are now given by 

£mn(%0i ^-l) 



exp 



(n + u ) (ni+ui) 
2m{ x H Xi} 



Li 



(3.11) 



(3.12) 



Note that this means that the £ mn (xQ,xi) are not periodic, being sections of £, but have 
the holonomy predictable from the character x defined in 3.8 above. The eigenvalues can 
be easily found and are given by 



A mn = ( 7— ) [l r l 2 ( n ° + u °) 2 ~ 2 ( n ° + u o)( n i + u i) T o + ( n i + u i) 2 ] 



\LqtJ 
( 2tt 

\L Tl 



(3.13) 



(n + ui) - r(m + u )|' 
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A computation of detA£ by the ^-function method establishes (3.6) where Tg(T 2 ,£) is 
given by (3.9). 

This completes our brief description of the continuum spectral data and we pass on 
to the matter of discretisation where we shall be able to reproduce exactly the continuum 
partition function — the inverse of the <9-torsion — in calculating the finite size effects on the 
lattice. 

To this end we now replace the r parallelogram by a discrete lattice which we then 
triangulate as follows: To retain complete generality over the geometry of the triangulation 
we construct a triangular lattice composed of similar triangles, pairs of which form paral- 
lelograms, The basic triangles have two sides of lengths ao and a\ with an angle 9 between 
them. The complete lattice forms our torus, T 2 , and consists of K G Ki sites and 2K Q Ki 
triangles, forming a parallelogram of sides L = K ao and L\ = K\a\. 

The resulting lattice is depicted in fig. 1 and we can use this figure to deduce that 
with this model for the torus a point z G C is identified with points of the form z + hiLq + 
nLi exp[i9] where m,ti6Z. Hence if z = (x, y) then 



(x, y) = (x + mil cos 9 + nLn, y + mL\ sin 9), m, n e Z 



(3.14) 



This geometric information also determines the holonomy of a section ip of the bundle C 
over T 2 which we record in the equation 

<p{x + mL 1 cos9 + nL ,y + mL 1 sm9) = e 27ri{rnUl+nuo) (p(x,y), m,neZ (3.15) 

We shall take a statistical mechanical standpoint and compute the energy of the lattice 
configuration rather than the, equally possible, approach of computing its action. To do 
this, we label the lattice sites by (ko,k\) = k = ko + (ki — 1)K with hi = 1, . . . , Ki, and 
the energy is then given by 

S C [T 2 , <p\ v] = \Y, v^*(k) (A(k, k') + m 2 5 k ,k') v>(k') ( 3 .1 6 ) 
kk' 

where ^fg = ciocii sin 9. 

The discrete eigensections of A are S noni (kQj k\) where 



£n m(ko, h) = exp 



K Ai 



(3.17) 



and the discrete version of the Laplacian A is a K Ki x KqKi symmetric matrix A(K 0: Ki) 
a general element of which we denote by A{(/c , &i), (k' Q7 k[)} so that we have 



A(tf ,tfi) = [A{{k M)Ak',l')}] KQK ^ KoKl 
with 1 < k , k' < Kq, and 1 < fei, k[ < K t 



(3.18) 
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Ki 

The partition function that interests us is 

n 

(n ,ni) 



and W = - InZ. 

We pass now to the lattice determinant in which we are interested namely 



(3.19) 



Using a nearest neighbour interaction all the non-zero entries of A(Ko, K\) can be deduced 
by its symmetry and by recording explicitly that 

A{(fe , fei), (k + 1, fei)} = -a = — \\ - 

sin 2 9 \ a o a o fl i / 

a r/7 7 \ /7 , -m r, 1 /I COS#\ 

A{ fe , fei , (/c , fei + 1)} = -/3= -—2- — 

sin 9 \ a i a o«i / 

A{(fe + 1, fei), (ko, h + 1)} = - 7 = (—) 

sin 6* \ a o a i / 

A{(ko, fei), (feo, fei)} = 2a, where cr = a + /3 + 7 

the remaining elements being zero. The lattice eigenvalues are then 

Anom = 2 (cr - acos(x no ) - /3cos(x ni ) - 7cos(x no - x ni )) 

, 2ttK + ^) . (3.20) 
where x ni = — , % = U, 1 



ZkqK! {T 2 , £, m) = / [n^'H = J[ { ^J^, } (3.21) 



det /A(k,k-) +m y k A /W"') ( 3 . 22 ) 

The computation of this determinant is a non-trivial task but one that we shall accomplish. 

Purely for convenience of the reader following the calculation we set K = K\ = K — all 
our results can be derived without making this specialisation and we shall quote formulae 
later (in the next section) without this restriction imposed. Now we want to compute 
ri(no m)(An ni + ^ 2 )/2 and it suits us to compute the logarithm of this quantity which we 
shall denote by lndet #-((A + m 2 )/2) for short, i.e. 

lndet K (^)= <A E 1, ,„(^^!) (3.23) 

It turns out that we can actually make sense of this expression and even evaluate its 
K — > oo continuum limit where we shall find a modular invariant expression. A key step 
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in achieving this is that we do one of the sums in 3.23 completely. We now outline how 
this latter step comes about. 

Let us define 5, by 5 = a + m 2 /2, so that using 3.20 this implies that 

Anom + m 2 = 2 (5 - a cos (x no ) - (3 cos (x ni ) - 7 cos (x no - x ni )) (3.24) 

Define the new quantities, b ni , z ni and 9 ni by writing 

, , r . , x ,, r .„ , (5 — (3 cos(x ni )) 

6 ni = (a + 7exp[-ix ni J) = |6 ni |exp[i0 ni J, z ni = — 

|"ni I yo.2o j 

=>• [5 -acos(x no ) - pcos(x ni ) - 7cos(x no - x ni )] = \b ni \(z ni - cos(x no + ni )) 
Then we factorise (z ni — cos(x no + 9 ni )) by observing that 



(z ni - cos(x no + 9 ni )) = ^ (1 - z- exp[i(x no + ni )])(l - z_ exp[-i(x no + ni )]) 
where z T = z ni =F yj - 1 

(3.26) 

Hence the sum over n in 3.23 requires us to compute 

(K-l) 



n ,ni=0 
(K-l) 

= In I b ni I (z ni - cos(:r no + 6 ni )) (3.27) 

n ,ni=0 

(K-l) ,, . (K-l) 

= * £ ln(-±i£-li) + £ ln|l- 2 _exp[z(* no +# ni )]| 2 



ni=0 no,ni=0 

But 

(K-l) (AT-l) 00 



ln(l - z- exp[i(x no + 6 ni )]) = - ^ exp[ir(x no + ni )] 

n =0 "0=0 r=l 

2^ ex Pb ra n ] 

r=l ' "0=0 (3.28) 



-E (z - e ^^" J) '" g^P""°'] 

1=1 



IK 



= ln(l-<) 
where v ni = z_ exp[z6> n J exp[27u-u /-^] 
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Now( c ) we note that we have entirely done one of the summations. It turns out that we 
can do still more: for the first term In (z + \b ni |/2) we can accomplish the remaining sum 
over rt\. We have 



In 



z+\b ni 
2 

where A 



ni 



But using the integral identity 

^ 2 doo 



n 



In 



1 A 



+ 



Al 



2 V 2 V 4 

26 — 2(3 cos(x ni ) 



ln(^ 2 -A 2 cos 2 H) = ln 



z + \Jz 2 - A 2 



(3.30) 



(3.31) 



we find that 



In 



z+\b ni 



^ 2 du 



In 



o 



7T 



A 2 

— - - |6 ni | cos {u / 



Now we have to do the summation and evaluate 

1 i^o Jo - 



(K-l) 
m=0 



In 



|6 n J 2 cos 2 (w) 



(3.32) 



(3.33) 



Next we reuse the two properties summarised in 3.26 and 3.28 to do the summation. But 
first we must rewrite 3.33 in a suitable factorised form: noting that A ni is quadratic in 

cos(x ni ) we get 



A 2 



\b ni \ 2 cos 2 (u) = [3 2 (cos 2 (x ni ) - 2qcos(x ni ) + 2a) 
= (3 2 (cos(x ni ) - x°)(cos(x ni ) - x 1 ) 



with 



x° = q 



+ vV - a 2 
x 1 = q — y 1 q 2 — a 2 



So 3.26 



A 2 



\b ni \ 2 cos 2 (u J )=(3 2 ^\l-x _\ 2 \l-x 1 _\ 2 



where i 



xl- = x 1 =F ^(x 1 ) 2 - 1 



( c ) In 3.28 above we use the fact that ex.p[2irir / K] are the K th roots of unity so that 

(Jf-i) 



^ ^ exp[27rmr/iC] = ^ 



K, if n = mod K 
otherwise 



(3.34) 



(3.29) 



r=0 



a similar algebraic trick was used in [11]. 
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If we now employ 3.28 we find that we can do the sum over n\ leaving us with a trivial 
two term sum originating in the roots of the quadratic. We get 



i \K\2 



2-KlUx . 

, % = o, 1 



where s l _ = x l _ exp 



K 



Thus we have accomplished our summation goal and have found for lndet^ ( A + m2 j 
the expression 



, J fA + m 2 \ . r^du, f(3 2 x° + x\\ * , , iN 



K\2 



i=0 
(K-l) 

+ E hi-<i 2 

ni=0 

(3.36) 

Actually the coefficient of K 2 above (the bulk coefficient), is easily checked to be given by 

7 7 1 4^) = y (2tt) 2 ^ - acos(^i) -/3cos(z/ 2 ) -7cos(^i - t/ 2 )}] 

(3.37) 

We can now begin to calculate what happens when K is large. Note that we shall 
now (temporarily) specialise to the case where m = since our answer in this case should 
be checkable against the B -torsion Tg(T 2 ,C) . However after our successful check against 
Tq{T 2 ,C) we shall give the result for m ^ also. We shall start with the two terms 
In |1 - (s l _) K | 2 , % = 0, 1 it turns out that 

lim K /" r/2 dU In If (s i ) K \ 2 - { °' for z = because \s°_ \ < 1 for u G [0, tt/2] 
k^oo J 7T _ \ 7^ 0, for z = 1 because |sL| < 1 for u G [0, tt/2] 

(3.38) 

We briefly describe the calculation of the surviving limit in 3.38 above. The entire 
contribution to the limit comes from the u = end of the integral since it is only at u> = 
that |si| attains the value 1. Taking this as our starting point we just need the (laboriously 
obtained) fact that, near u = 0, 

| T |2 

xL = l-— uj + --- (3.39) 
which is relevant since si = x\_ exp[2%iui/K]. This found we write 



I 



K,2 —\n\l-(s 1 _) K \ 2 = r—\n\l-{s 1 ) K \ 2 + f /2 — lnll-^)^! 2 , e>0 (3.40) 
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The second integral dies exponentially asX^oo but the first integral gives, for small e, 



1 \K\2 



K 



f-Hl-C* 1 ) 
Jo n 

r—\n\l-{s 1 ) K \ 2 
Jo 71 



1=1 



t^S~^ f e du exp[27rmi/] / |r| 2 

h'° ~ — ' — ' ~ 



2 - » ! ' 
■UJ 



^ exp[27riui/] / -ri \ / |r| 2 \ 



2 \ ffl+r 



From which we readily compute that 

^ 2 dw 



(3. 



o 
1) 



/" doj i 2n v-^ COs(27rUiZ) 
hm K — \nl-(s 1 _) K2 = , 2 

K^oo 7 7T 7T|r| 2 ^ / 2 



27TTi 



(3.42) 



Having taken care of the second term in 3.36 we move on. 

The first term in 3.36 already has its independence displayed as an explicit factor of 
K 2 and so the only other term to consider is the last term , namely ^i^=o hi |1 — v n 1 1 2 - 
We note that, apart from an explicit factor of K, the dependence on K in the sum is via 
the variable 

2-7r(ni + u\) 



K 



(3.43) 



and for K — > oo we distinguish the two cases 



K — > oo 



ni/K — > 0, case 
ni/K y4 0, case 1 



(3.44) 



Since the procedure for both cases is lengthy but reasonably straightforward we describe 
it fairly briefly: first we deal with case and then summarise the differences relevant for 
case 1. 
Case 0: 

We consider first the term In 1 1 — \ 2 and simply quote the fact that, after slightly 
tedious algebra, one finds that, with A = \J (a[3 + ]fy + 7a) / (7 + a), we have 



exp 



27rA(ni + u\) 



K 



1 - 



Z7 2ir(ni + ui) 



(a + 7) 



exp 



-27TT1 

K\t\ 2 



ir 27r(ni + ux) 
Irj 2 K ,C ' X1) 



K 

K 



exp 



K 



(3.45) 
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Hence we now obtain 



lim v„ = exp 



-2m \ (m + Mi) i— To - u o 



exp 



. (ni +ui) 
—2m { u 



(3.46) 



This expression 3.46 will take care of the K — > oo limit of the term In 1 1 — | . We must 
now round off the treatment of the K — > oo limit by describing what happens for case 1. 
Case 1: 

Recall that case 1 means that ni/K as K — > oo. We can achieve this by writing 



m = (K-l), l = l,...,(K-l) 



(3.47) 



But since v ni only depends on n\ via the functions cos(x ni ) and sin(x ni ) it is enough to 
note that with m h- > (if — /) we obtain 



cos(27r(ni + u±)/K) = cos(x ni ) i— > cos(2-7r(/ — u±)/K) 
sin(27r(ni + Ui)/K) = sin(x ni ) h- > — sin(27r(/ — ui)/K) 



(3.48) 



The discussions of case and case 1 are now complete and this allows us to deduce that 



lim 



t>ni = ex P 



. (/ — Mi) 

2m { ± — — ^ + u 



We now have the entire K — > oo limit of lndet (A/2) and it is that, as K — > oo 
'A 



(3.49) 



In det k 



o /"^ db'ldb'2 , r r / \ ^ / \ / , 

K I — — In [{cr — a cos(z/i) — pcos(z/2) — 7 cos(z/i — z/ 2 )}J 



(2ny 



2ltTi i 2 



+ ^ In 1- exp -27rz | + e n - w 



where e n = 



1, n > 
-1, n< 



Hence the finite part of lndet kA/;/2 is the function F(it , «i, r) given by 



(3.50) 



F(uo,u lt T) = — ^ (u\ -m + ^ J + ]T In 



1 — exp 



-2m { — + e n ^ — - u 



(3.51) 

Now to agree with the continuum calculation F(uo, u±, r) should give the <9-torsion via the 
equation 

\nT B (T 2 ,jC) = F(u ,u 1 ,r) (3.52) 
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But T B {T 2 ,C) was given in 3.9 from which we obtain 



In T B {T 2 ,C) = In 



ex.p[niu r 



2 0i (ui - tu ,t) 



t]{t) 



2 



= -2ttti ( ul - u + - j + ^ ln|l - exp [27rz {|n|r + e n («i - w t)}]| 2 , 

(3.53) 

(The last line of 3.53 follows from the series representations of rj(r) and 9\{u\ — tuq,t).) 
However the apparent disagreement between 3.52 and 3.53 is illusory these two expressions 
are both modular invariants and if we perform the summations leading to 3.52 in the 
opposite order then we get precisely 3.53. Hence we do indeed find that the continuum limit 
of our lattice model has reproduced the <9-torsion Tq{T 2 ,C) in a fashion quite analogous 
to the one dimensional example of the previous section. 

It is useful at this point to spell out some of the details of the modular properties of 
our expressions in a short discussion. We recall some basic facts. If M denotes an element 
of the modular group then M acts on the modular parameter r according to the well known 
rule r i— > (ar + b)/(cT + d) with integral a,b,c and d satisfying ad — be = 1. Hence we usually 
write something like 



M=(^ c e 5L(2,Z)/Z 2 , a,6,c,deZ, ad-cb=l 



(3.54) 



However when such transformations M act on the torus they act on cycles and thus on the 
generators of the fundamental group 7Ti(T 2 ); or equivalently on the holonomy phases Uq 
and u\ which specify the bundle C. Hence we must work out this action on the phases. But 
the modular group SL(2, Z)/Z 2 is generated by the celebrated S and T elements which are 

S:t\ — ► --, T:n — >t + 1 (3.55) 

T 

Therefore it is sufficient to give the action of S and T on the holonomy phases and this is 

Under S : uq i — >■ Uq, U\ i — ► Uq + U\ 

(3.56) 

Under T : uq i — >■ u\, U\ i — >■ — Uq 

In addition S and T are isometries of the flat torus so they leave the volume V invariant. 

With this information on S and T it is possible to check directly the modular invariance 
of 3.52 and 3.53. For completeness we give the corresponding information for a general 
element M of the modular group. It is that, regarding L and L\ as complex numbers in 
an obvious way, we have 

Under M : uq i — > cu\ + duo, u\ i — ► au\ + buo 

and t i — > aT — (3.57) 
ct + d 

with M G SL(2,Z)/Z 2 



17 



Modular invariance, lattice field theories and finite size corrections 



Nash and O' Connor 



For the 9-torsion this modular invariance is the assertion that Tg(T 2 , C) only depends 
on the conformal structure on T 2 i.e. on its complex structure. In fact as observed in 
[10] modular invariance in this case can be shown to follow from the classical functional 
equations obeyed by 9\ and r\. 

In fact in this work the modular invariance of the continuum theory is guaranteed by, 
and has its origin in, the independence of the limit defining F(w , , ui,r) on the order in 
which its double sums are carried out. 

Finally we reinstate the mass and give the corresponding formulae when m/0. We 
find that 

lndetx ^ A ^ m ^ _> k 2 J ^|y^ In [{5 - a cos(z/i) -/3cos(z/ 2 ) -7cos(z/i - v 2 )}] 



tyti m z V . , 
— c[uq, ) + > In 

D Tl * — ' 

n= — 00 



1 — e 



(3.58) 

where V = L Li sin 9 and the function c(u, x) that appears in 3.58 above is defined by the 
equation 

f 00 * In 1 - 2 = - C ^^ (3.59) 

One can verify that 3.58 reduces to 3.53 when m = — for example this would require that 
c(u, 0) = 12(u 2 — u + 1/6) which is true. 

Of course, for this case where m^Owe can use a zeta function method to give the 
continuum partition function. Repeated use of (2.10) enables one to calculate (' d * dc+m2 (0) 
which is the required object. The result is that we obtain Wp in agreement with 3.58 
above; we also find a bulk term of the form — ^^-(ln[(m/27r//) 2 ] — 1) where ji is an arbitrary 
undetermined scale. 

We would like to emphasize that, in the continuum or scaling limit, we have found 
the exact finite size corrections to lnZ(T 2 , £, m 2 ) and these are modular invariant whether 
m = or not. But when m is zero we have a critical phase and then the finite size corrections 
to Z exhibit holomorphic factorisation and are given by elliptic theta functions, i.e. 

exp[-W F ] = TT 

r • 2 ^(m-TUQ.T) (3-60) 

with J- = exp 7rzM n r — 

7/(r) 

Concerning holomorphic factorisation note that if the phase factor exp^m 2 , 7 "] were 
absent from JF above then Z would also exhibit holomorphic factorisation in the Picard 
variety label w = u\— tuq. Actually the presence of the phase exp^i-u^r] is a consequence 
of the central charge c being non-zero — lack of holomorphic factorisation in this sense is 
the existence of a holomorphic anomaly [12] of the Be determinant bundle over the Picard 
variety. If we observe that 

TV 

d w dw \nTg(w, w, t) = 4ir5(w, w) (3.61) 

Tl 
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we see that the torsion is proportional to the Greens function, where the lowest mode (the 
zero mode) has been projected out, of the Laplacian d w du, on the Piccard variety. It is the 
absence of this lowest mode that is the barrier to holomorphic factorization. 

In the next section we turn to the bulk term and the cylinder limit L\ — > oo. 

§ 4. Conformal properties and effects due to finite size and finite lattice spacing. 

In this section we work with Kq, K\ and the mass m at perfectly general values. 
We recall that W and Z are related by 



W = -InZ 



(4.1) 



If we now perform the lattice sums, doing first the sum over no, followed by that over n\ 
we obtain the result that 

W = K Q K{W B + W F (4.2) 

where 



W B = 



w dv\dv2 



In 



n 



{5 — acos(^i) — /3cos(z/ 2 ) — 7cos(z7 — z/ 2 )} 



and 



W F =K i y ^ln|i_( a *_)^o|2 + y ln |i 



(4.3) 



(4.4) 



We note that Wb gives the free energy per lattice site in the thermodynamic limit, i.e. we 
have 

W 

lim — — = W B (4.5) 

Ko^^oa K K 1 

This in turn means that Wf gives the complete finite size corrections to the bulk lattice 
behaviour. We now examine in more generality and more detail the limits studied in the 
preceding section. We now no longer set Kq = K\ = K we simply take Kq,Ki — > oo 
while keeping fixed their ratio k = Kq/Ki, we also keep fixed the quantities m 2 , 9 and 
Li = Kocii, we shall refer to this as the scaling limit. We still find the same result for the 
finite size corrections Wf namely 

2 



lim Wf 

scaling 



2t r 00 

•— c(u , )+ V ln 



t -2-KT 1 ^J(n+u ) 2 + ^^+27Ti{u 1 -T (n+u )} 



while for the bulk term we find the mass dependence 

Vm 2 



6) 



lim K K 1 W B = K K 1 K B - , 

scaling 47T 



Vm 2 m 2 V 

MKoKi] - 2p} - — (M^zd - 1) + • 



47T 



4tt 2 



where A B = W B \m=o, 



and p = p(a, (3, 7) = / dv 

Jo 



V = L Li sin 9 
1 



sin v\f \ + go 2 sin 2 v 



-ln[^(/3 + 7 )] 



(4.7) 



19 



Modular invariance, lattice field theories and finite size corrections 



Nash and O' Connor 



Our results can be applied to the Ising model on a torus — the special case correspond- 
ing, in our language, to tq = was studied in [13]. This arises because of the equivalence 
of the Ising model and a dimer model on a decorated lattice cf. [14,15]. To see what 
happens we denote Wf by Wf{uq,u\) in order to display its phase dependence and use 
[16] for a, (3, 7 and 5. Let us do this: denoting nearest neighbour spin-spin couplings by 
Ji, J<i and J3 then 



a = sinh(^M, (3 = sinh(^L), 7 = s inh(^;) 

Kb-L Kb-L kB-L 

x u/ \ , / 2J2 , , 2J 3 . 2Ji . 2J 2 . 2J 3 

5 = ^kjf) COSh( ^T ) C ° Sh( ^T ) + ^kjf) ^kjf) 



(4.8) 



and the Ising partition function is given by 

This result is for ferromagnetic couplings, with + referring to T < T c and — to T > T c . If 
we take the the scaling limit Wp(uo, u\) 1— > Tf(uq, u\) and we obtain a modular invariant 
expression for the finite size contributions in the scaling limit for this model also. Our 
results (4.9) incorporates the complete lattice and finite size corrections for the Ising model 
on a triangular lattice. With a similar equivalence our results can easily be translated to 
give the general result for other models. 

We point out that the expression of Z Isirig as the four term sum 4.9 above constructs 
a modular invariant function of r and m 2 V; this construction makes it easy to understand 
that it is the summing over the ttj which form an orbit of the SL(2, Z)/Z2-action on the 
space of iii's that guarantees the invariance. Mathematically we are working with the 
action of SX(2,Z)/Z2 on the space of flat bundles C — the so called Picard variety — and 
many other orbits exist, this allows us to construct many additional (phase independent) 
modular invariant partition functions such as those of the other conformally invariant field 
theories cf. [17] and references therein. Rational conformal theories will be obtained when 
the Ui orbits contain elements corresponding to roots of unity. 

Next we would like to consider an interesting geometric limit of the model where 
L\ — > 00: the cylinder limit. This will enable us to access the central charge of the model. 
For large L\ the quantity Wf/V tends to the finite value^ Y yhnder where 

^cylinder = _*^ UOfm 2 L 2 ) (410) 

We refer to c(w ,m 2 -^o) as the 'cylinder charge', the central charge c of the model is 
obtained by setting m = uq = and we see that this gives us 

c = 2 (4.11) 



(d) 



We could equally take Lq large if we interchange (uq, Lq) with (ui, L\). 
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Because the cylinder charge determines the central charge at appropriate values of its 
variables it might appear reasonable to expect that the cylinder charge ^cylinder be equal 
to the Zamolodchikov c-function [18] — this is not so. To see this note that the c-function 
must be monotonia in m and the cylinder charge is not cf. fig. 2 where we plot c(u, x) versus 
x for various u. In fig. 2 we see that c(u, x) exhibits a universal crossover phenomenon as 
x — > oo with c(u, x) having the property that lim^oo c(u, x) — > whatever the value of u. 

The cylinder charge for the Ising model, can be obtained from 4.9 with L\ — > oo, and 
it works out to be 7/^^ er where 

^i:t r = -^^x) (4.i2) 

This means that comparison with 4.10 above determines that the cylinder charge for the 
Ising model is 

( 4 - 13 ) 

Now if we set x = we should get the Ising central charge and doing this we find c = ^ 
(cf. fig. 2) the conventional value [19,20]. 

A further rich feature of this model is that is possesses non-commuting limits namely 
the limits uq,u\ — > and m — > 0. To establish this we can expand the finite size Wf 
contribution to 3.58 for small U{ and m obtaining 

W F = In [(2n) 2 \ Ul - tu \ 2 + r im 2 V] + 2 In |??(t)| 2 + • • • (4.14) 

But the RHS of 4.14 can tend to the two distinct (logarithmically singular) expressions 

ln|ui -rw | 2 + 21n|?7(r)| 2 
and ln[rimV] + 21n|?7(r)| 2 

depending on the order in which the limits are taken. Nevertheless both limits and indeed 
4.14 itself are modular invariant. 

This model is also ideal for studying the Kosterlitz-Thouless phase [21]. One can 
add on a |</?| 4 interaction and in this way one can study the critical point of the XY 
model approaching it from the disordered phase. We have holonomy conditions (3.15) 
that describe the essential features of a vortex phase: To make this more transparent one 
should perform the cylinder limit of 3.15 and then map the cylinder to the plane using the 
conformal map 

y + t— — =lnz (4. lb) 

This gives a vortex at the origin whose charge is given by uq and whose presence is detected 
by the existence of an Aharanov-Bohm effect for transport around the origin. 

§ 5. Conclusion 

We have considered lattice models and their continuum limits in one and two dimensions. 
Our one dimensional model is that of a gauge theory for a flat bundle L over the circle and 
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it is non-trivial its partition function being a power of the Ray-Singer torsion T(S 1 , C) of C. 
We find that the combinatorial and continuum partition functions agree precisely and are 
equal to T(S , This result emphasizes the topological nature of T(S , which is 

of course expressible purely combinatorially, as it was originally [22], or analytically as it 
was later [23]. 

In summary, the finite size corrections to the free energy are modular invariant. This 
conclusion extends to the entire scaling neighborhood of the critical phase. We use our 
results to give expressions for the complete lattice and finite size corrections for the two 
dimensional Ising model on a triangular lattice via its equivalence to a sum over Pfaffians. 
Modular invariance also extends to models in more than two dimensions when the geometry 
giving rise to finite size effects contains a flat torus. For a three dimensional cylindrical 
geometry with toroidal cross-section the result can be obtained from Wp by replacing m 2 
with m 2 + q 2 and integrating the resulting expression over q. One can understand the origin 
of modular invariance in general as the residual freedom to reparametrise coordinates, in 
the continuum limit, while retaining flat toroidal geometry. 

In the two dimensional case the limiting finite size corrections at the critical phase are 
expressible in terms of classical elliptic functions. Infinitesimally small values of the phases 
Ui lead to logarithmically divergent contributions to the free energy. This means that the 
free energy needed to create a vortex becomes infinite for an infinitely large lattice. In 
general the model has a surprisingly rich structure of non-commuting limits. For example 
the limits of approaching the critical phase (m — » 0) and that of sending the Ui to zero do 
not commute. 
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The triangulated torus 




Fig. 1. 

The 'cylinder charge' function c(w, x) for various u 



c (u, x) 




Fig. 2. 
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